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I. INTRODUCTION 



Putting an electron into a crystalline lattice, it will interact with the dynamical deformations of that latter. Depend- 
ing on the type of material we are dealing with — such as ionic polar crystals, covalent materials or simple metals — the 
form of the electron-lattice coupling as well as the relevant lattice modes being involved in this coupling will be differ- 
ent. The mutual interaction between the electron and the lattice deformations, which in general is of a highly complex 
dynamical nature, results in composite entities of electrons surrounded by clouds of virtual phonons which, in a dynam- 
ical way, correlate the position of the electron and the associated local lattice deformation. Such entities are referred 
to as polarons, which can be either spatially quite extended "large polaron", or rather constraint in space "small 
polarons" . Large polarons generally are itinerant entities, while small polarons in real materials, have a tendency to 
self-trap themselves in form of localized states. By far the most spectacular manifestation of such electron-lattice 
interactions in the weak coupling large polaron limit is the classical phonon mediated superconductivity. 

The earliest indications for polarons came from F centers — electrons trapped in negative ion vacancy positions in 
alkali halides — which led Landau to the concept of localized strong coupling polarons 7 . Subsequently, the field was 
developed mainly in connection with scenarios where the materials could be characterized by a continuous elastic 
medium and specifically for: 

(i) ionic crystals 7 , where positively and negatively charged ions oscillate out of phase around their equilibrium posi- 
tions and thus give rise to large electric polarization fields in form of longitudinal optical modes, which then strongly 
couple back onto the conduction band electrons, 

(ii) covalent materials 7 , where local dilations of the material couple the electrons to corresponding deformation po- 
tentials and 

(iii) simple metals 7 , characterized by high frequency plasma modes which, because of strong charge density fluctua- 
tions, renormalize into acoustic modes. 

With the arrival of new materials since the late sixties, the necessity to consider the microscopic lattice structure 
for the polaron formation became more and more crucial and initiated the theoretical work on the transition between 
large and small polarons and the question of polaron localization. 

After reviewing the different kinds of electron-lattice coupling leading to polaron formation we shall discuss the 
fundamental issues of the cross-over between large and small polarons, the question of continuous versus discontinuous 
transition and the difference of the polaron self-trapping with respect to localization in systems with attractive 
interaction potentials. We then present examples of a few decisive experiments on polarons which, early on, shed 
some light on their dynamical formation and disintegration. 

II. THE FROHLICH LARGE POLARON 



A. Polarons in ionic crystals 



Let us to begin with consider the case of an ionic crystal where longitudinal optical phonons strongly couple to the 
electrons. A single electron in such a dielectricum induces via its charge an electric displacement field 

D(r, r c i) = -grad— , (1) 

l r r el| 

at a spatial coordinate r with r c i denoting the position of the electron. This displacement field couples to the dynamics 
of the lattice which can be described in terms of a polarization field P(r) of the dielectric medium and which (via 
Poisson's equation div D(r) = 47re(5(r — r c i)) can be expressed in terms of a polarization potential $(r) such as: 

P(r) = -Igrad $(r), $(r ol ) = -- [ d 3 r D(r, r cl ) • P(r) . (2) 
47r e ./ 
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We assume for simplicity the absence of any shear and vorticity in the elastic medium forming the dielectricum. 

The standard way to describe the physics of polarons is to introduce measurable material quantities such as the 
frequency dependent dielectric constant s(lo), defined by D = e(w)E, in order to account for the dynamics of the 
dielectric medium via the frequency dependence of the polarization field, related to E via 47rP(r) = D(r) — E(r). 
D(r) being exclusively determined by the point charge is independent on e{uo). In a crude sense one can consider 
the response of the dielectric medium as being given by a superposition of two contributions: (i) a high frequency 
contribution in the ultraviolet regime £oo, arising from the electron clouds oscillating around the ionic positions and 
(ii) a low frequency contribution e m the infrared regime, arising from the oscillations of the positive and negative 
ions against each other. This separates the polarization field P(r) into two contributions P(r) = Po(r) + P OG (r) which 
in the low and high frequency limits are determined by: 

P„(r)+P„(r) = i-(l-£)D(r), (3) 



considering the fact that the low frequency contribution of the polarization field remains unaffected by high frequency 
perturbations. The determinant contribution to the polaron dynamics, arising from the low frequency polarization 
field is hence given by 

Before entering into a detailed discussion of the intricate nature of the dynamics of this problem let us consider 
certain limiting cases and start with the picture of an electron inside a dielectric continuous medium being constraint 
to a finite small volume in a sphere of a certain radius R\ to be determined. 

We consider an electron, moving inside such a sphere of radius Ri with a velocity v, to be fast compared to 
the characteristic time of the atomic oscillations (2tt/loq) (given by the longitudinal optical lattice modes). The 
polarization field induced in the medium by the motion of the charge carrier can then be considered as static for 
distances much greater than 2ttv/ul) and being described by a static Coulomb potential. Inside this sphere the test 
charge can be considered as being uniformly distributed and hence the potential as being constant. We thus have: 

E pot = -e 2 /iR 1 (r < J2i), E pot = -e 2 /er (r > J2i) . (6) 

Determining the size of the polaron in a purely static fashion, considering the electron kinetic energy inside this 
sphere as being given by Skin — {2-kK) 2 /2mR\) 1 we find by minimizing the total energy with respect to R\\ 

R 1 = {2irh) 2 e/me 2 . (7) 

Determining the size of the polaron radius in a dynamical fashion we require that the characteristic wavelength 
2irh/mv of the electron must be smaller than 2-KvjujQ. This implies that the typical radius Ri of the trapping potential 
is determined by the distance for which h/mv ~ v/ujq, or in other words: 

R 2 = 2ir(h/mujo)^ . (8) 

^From these simple arguments it follows that the interaction between the displacement field and the polarization 
field, arising from the dynamics of the lattice, (the polar coupling between the electron and the longitudinal optical 
phonons)) is of the order of E pot — — e 2 /ei? 2 - 

The potential energies for the two limiting cases where the dielectric medium is considered (i) as static and (ii) as 
dynamic, are given by 

_ t/l _ {U 2 V _ 1 2 

where apv is a dimensionless coupling constant introduced by Frohlich 



=2 
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The value for a^r can vary over one decade in the interval ~ [0.1, 5] depending on the material. Large values of a favor 
the static picture, while smaller ones favor the dynamical one. For most standard materials (polar and ionic crystals) 
the typical phonon frequency is of the order of 10 -14 seconds and the electrons can be assumed to be essentially 
free electrons. This results in a value for i?2 which is much bigger than the lattice constant and hence justifies 
the continuum approach described above permitting to treat the dielectric constant as wave vector independent. 
Within such an approach Frohlich first formulated this polaron problem in a field theoretical form and proposed a 
corresponding Hamiltonian for it, designed to describe the so called large polarons. 

Let us now briefly sketch the derivation of the corresponding Hamiltonian 7 . To within a first approximation the 
dynamics of the polarization field can be described by harmonic oscillators driven by the electric displacement fields. 
Dealing with the polaron problem on the basis of the continuum model of a dielectricum one is restricted to consider 
waves of the polarization field with wave length much bigger than the inter-atomic distance and moreover to limit 
oneself in this problem to frequencies well below the optical excitations which would result from the deformations of 
the ions rather than from their motion. Under those conditions one can limit the discussion of the dynamics of the 
polarization field to the infrared longitudinal component of it, being driven by a source term due to its coupling with 
the electric displacement field. Using a simple Ansatz in form of a harmonic motion of the polarization field one has 

,l ~ A P (r) = iD(r,r cl ). (11) 



dt 2 1 -<y- u ^' 7 J 

2 

The coefficient i = = of D(r, r i) on the rhs of this equation is determined from its static limit, following eqs. 10 
. To this equation of motion for the polarization field one has to add the term which describes the motion of the 
electron and which is controlled by the interaction energy e<I>(r), given in eq. J5J), i.e., 

d 2 

m—r c \ = -e grad $(r c i) . (12) 

Considering 7 Po(r), mv c \ and 7Pp(r), p e i = mr e \ as a set of generalized coordinates and conjugate momenta 
{Q, dLdQ} the corresponding Hamiltonian is given by 

H = lyj d 3 r [Pg(r) + wgP§(r)] - j d 3 rD(r, r ol )P (r) + \mvl . (13) 

In order to obtain a quantum field theoretical formulation of this problem one introduces habitually the vector field 
representation for the polarization fields 



qA 

with the phonon annihilation (creation) operators a q ^ with [ciq\, a^ry] = fiq.q'Sx.y ■ 

For ionic polar crystals the dominant contribution of the polarization field comes from longitudinal optical phonons. 
We shall in the following restrict ourselves to those modes only. Introducing the quantized form of the electron 
momentum p i = with [p e i,r e i] = — ih, we rewrite the Hamiltonian, eq. I|13l) in terms of the phonon creation 

and annihilation operators and the electron coordinates and conjugate momenta as 

ff =f| + ^°E (< a i + I) + E V^ia+e-^ - a q e^) , (15) 

where the effective electron lattice coupling constant is 

U q opt = iV4^{h/2Lu m)ihuJo(l/q) . (16) 

Upon introducing dimensionless electron coordinates as well as phonon wave vectors r c i = (2muj /h)^r c i and q = 
(2mLOo j K)~ "5 q, H/Hluq turns out to be exclusively controlled by the unique parameter aFr- 

Replacing in the expression for the Hamiltonian, eq. (|15|) . the Fourier transform e lq r of the electron density p(r) = 
^2 r i 5(r — r e i) by its second quantization form ^2 ka c£ +qcr Ckcr we finally obtain what is generally called the Frohlich 
Hamiltonian 

H = E £ * C i« C ^ + E ^0«flq + \) + E ^'i+qA^S - <»±q) • ( 17 ) 
kcr q qk(T 

£k = D — i7k, 7k = ~X) 1 5 ex P(*^ : ' denotes the bare electron dispersion (with S being the lattice vectors linking 
nearest neighboring sites) and D = zt is the band half- width and z the coordination number. 
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B. Polarons in non-polar covalent materials 

A structurally very similar Hamiltonian is obtained for covalent materials when the density of charge carriers is very 
low and screening effects can be neglected. The electron lattice interaction then can be derived within the so-called 
deformation potential method 7 where it is assumed that the electron dispersion of the rigid lattice e£ gets modified 
due to some elastic strain which, in the simplest case, amounts to a dilation. This then results in a shift of the electron 
dispersion = e£ + CA with 

a « - «f - « vsf £ Tjfe^*" - < e "-"-' ■ < 18) 

The dilation, being exclusively related to longitudinal acoustic phonons (a q , a^) with a phonon frequency w q and M 
denoting the ion mass, the constant C can be estimated form pressure measurements and is typically of the order of 
10 eV. The effective electron-lattice interaction in that case is then similar to that derived above for ionic crystals 
and is described by the Hamiltonian eq. El upon replacing V^ pt by 

V? = Q . (19) 

C. Polarons in metals 

A yet very different approach is required to treat the electron-lattice coupling in metals. The electrons at the 
spatial coordinate x experience a pseudo potential V(x— Rj) exerted on them by the positively charged ions situated 
at lattice sites R;. The major effect giving rise to the coupling of the electrons to the lattice now comes from the 
dynamical displacements = R; — R° of the ions from their equilibrium positions R^ . The dynamics of the lattice 
arises primarily from the ion-ion interaction V(Ri — Rj), which can being assumed in form of a Coulomb potential. 
But since those ions are imbeded in a kind of electron jellium, they sense the charge density fluctuations of the 
electrons via the electron-lattice interaction. Considering plane waves i/j(x) f° r the electrons one has the following 
Hamiltonian to consider 

H = jdxV,+ (x) -g+^x-R,) ^ (x )+£|jL + i£y(R i -R J .) 

Li J i i,j 



+ i / dxdx'V> + (x)^+(x')^^^(x>(x), (20) 



where Pf /2M denotes the kinetic energy of the ions at site i. Using the standard expansion of the lattice displacements 
Ui in terms of the phonon operators 

qA 



* = ^E^fe) (21) 



(e q ^ denotes the polarization vectors and f2 q A the eigen-frequencies of the various phonon branches) one can write 
the above Hamiltonian in the form 



kcr qA q 

+ ECA V q («qA+«i qA ), (22) 
qA 

where pq = ^ Sko- c k+qo- Cko ' 1S tne charge density fluctuation operator for the electrons and V^ * = 
— i(N/2MhQ,q\)^ (q • eq\)V(q). The relevant phonon modes arising from the ion-ion interaction are the practi- 
cally dispersion- less ionic plasma modes with frequencies f2 q A- However, these modes being coupled to the electron 
charge fluctuations, this results in a strong renormalization leading to: (i) a longitudinal acoustic branch with a 
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correspondingly phonon frequency w^ c =~ q^Jm/3MvF (vf denoting the Fermi velocity) and (ii) to a corresponding 
dressed electron-phonon coupling given by ? 



• N < ,A(q)VXq) 

v ^ = -^^u^-^^z7- (23) 

e(q) denotes a wave vector dependent dielectic constant, A(q) the Coulomb screened renormalization of the bare ion 
potential and Z c the spectral weight renormalization constant of the electron quasi-particle spectrum. 

In spite of the diversity of the physical systems considered above, the various Hamiltonians describing them are of 
rather general form. The main message which they contain is that, because of the electron-lattice interaction Vq, the 
electrons will be accompanied by a lattice deformation. In the weak coupling limit, this is restricted to simply a single 
phonon accompanying the electron. Within a lowest order perturbative approach this leads to states of the form 

<i°> + i:v t %xi°) ( ^ + ^ q - fk ) ' < 24 » 

where |0) denotes the vacuum state for the electrons as well as of the phonons. This expression clearly indicates 
that electrons carry with them an electrical polarization field or, in other terms, the electrons are accompanied by a 
phonon. This feature reflects the fact that the total wave-vector 

K = kc £r CkCT + J2 ^qA ^ ' ( 25 ) 

kcr qA 

is a conserved quantity. The most spectacular consequence of this lies in the formation of Cooper pairs 7 , as evi- 
denced by the isotope effect of their binding energy 7 and the ultimately resulting phonon mediated BCS theory of 
superconductivity 7 . Electron pairing then can be understood by a process in which the passage of a first electron 
polarizes the lattice and where subsequently a second electron reabsorbes that polarization. The effective interaction 
Hamiltonian for this pairing is given by second order perturbation theory: 

ffel-e. = £ |V q | 2 C+ +qT C+_ qi C kU C kT _ _^ _ . (26) 

kk'q ^ ^ 

This effective electron-electron interaction shows that within a small region around the Fermi surface i.e., |ek — £k± q | < 
w q this interaction is attractive and therefore leads to an instability of the Fermi surface resulting in a superconducting 
ground state via Cooper pair formation in k-space. 

This present discussion of the continuum approach to the coupling between the electrons and the lattice-deformations 
has shown interaction Hamiltonians which, depending on the underlying materials, show either a coupling between 
the phonon coordinates (a q + a+ q )/\/2 and the electron charge density p q = ^ X)k c k+s c k+a (metals) or between 

the conjugate phonon momenta (a q — a+ q )/i\/2 and the charge current density p q = qc k+a c k+a (polar and 

covalent materials). Formally these interaction terms can be written in a unifying way by rotating the phonon 
coordinate into the phonon momenta and vice versa by a suitable unitary transformation U = exp(— i^a+a q ). 



III. THE HOLSTEIN SMALL POLARON 



The systems discussed in the previous section are characterized by long range electron-lattice coupling which show 
up in form of (i) a moderate mass renormalization of charge carriers in their band states, (ii) an equally moderate 
reduction in their mobility due to the scattering of the electrons on the lattice vibrations and (iii) the emergence of 
phonon sidebands in optical absorption spectra. 

The fundamental theoretical question which posed itself in the context of polaron physics ever since Landau 7 
proposed self-trapped localized polarons was to establish whether in such systems one would have a local lattice 
instability upon increasing the coupling constant apr, passing from large mobile Frohlich polarons for weak coupling 
to localized polarons when the coupling strength exceeds a certain critical value. This question could not be addressed 
within the continuum approach since it requires a physics which is related to the dynamics of the local lattice 
deformations on the scale of the unit cell. On the experimental side, more and more new materials where synthesized 
in the mean time which clearly showed polaronic effects on such short length scales. It was for these reasons that a 
scenario was introduced which could describe such lattice polarons. 
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The generic model to capture such a situation, generally referred to as the Holstein molecular crystal model 7 treats 
the problem consequently in real space. Its corresponding Hamiltonian is given by 

H = Dj2n ia -tJ2(. c t c j* + h - c -)- X J2 UiUi 

+ ^y(»i+*i)+E [/ ' l it™il- (27) 



where m a = cf a c ia (ni = rii a ) denotes the density of charge carriers having spin a at molecular sites i. The 
electrons are assumed to be coupled to the intra-molecular deformations u\ via charge density fluctuations and the 
coupling constant is denoted by A. The dynamics of the lattice is treated purely locally with Einstein oscillators 
describing the intra-molecular oscillations with frequency luq. M denotes the mass of the atoms making up the 
diatomic molecular units. The additional Hubbard U intra-molecular repulsion is introduced sometimes in order to 
account for possible correlation effects in conjunction with the purely polaronic features. We shall not consider the 
effect of this term in this present discussion, but several specific lectures in this school will be devoted to it. 

This model is capable of describing the self-trapping of charge carriers which arises from a competition between the 
energy gain coming from the itinerancy of the electrons and that coming from the potential energy due to the induced 
local deformations of the molecular units. In order to see that let us rewrite this Hamiltonian in a form which makes 
more explicit such a self-trapped localized picture: 



i i=£jo i 



x - M 



Xrii \ 2 2 f Ani 



M 

T 



2uj- 



(28) 



The major features which evolve out of such a representation are: 



The energy of the electrons is lowered by an amount e p = X 2 /2Muiq which corresponds to the ionization energy 
of the polaron 



• A lattice induced intra-molecular attraction of strength 2e p between the electrons on different sites of the 
diatomic molecule units which can partially or totally compensate their intra-molecular Coulomb repulsion U. 

• The intra-molecular distance of the Einstein oscillators are shifted by a time dependent quantity u? —Xni/MujQ 
which follows the time evolution of the charge redistribution on site i. The frequency of the oscillators is modified 
as a consequence. 

Since in this school we shall primarily be concerned with systems which can be described by the Holstein model, 
let us now focus in more detail on the basic physics inherent in this scenario. To begin with, we consider this problem 
in a semi-classical fashion and in the adiabatic approach. 

In the limit of small coupling e v <C D, the form of the Hamiltonian, given in eq. (|28|l is reminiscent of one which 
describes itinerant electrons in a static potential, given by l/ wc ({wi}) = —X(ui)cf r c ia dfd i with (m) = X/Muj 2 and 
d^'d i representing the density of some fictitious particles. Solving the eigenvalue problem for the ground state energy 
E Q , i.e., 

A 2 1 

M^ ^^~E~o =1 (29) 

predicts a splitting off of the ground state energy from the bottom of the free itinerant band (E < 0), resulting in a 
localization of the charge carriers. 

In the limit of strong coupling, for large values for A such that e p > D, the adiabatic approach tells us to ignore the 
electron itinerancy and assume the electron to be fixed at a particular lattice site. This then results in an adiabatic 

potential for the electron given by V sc ({ui}) = —e p rii + ^-lu 2 (ui — jf^j > which follows directly from the form of 

the Hamiltonian given in eq. I|28|l . If one requires that this potential is deep enough to bind the electron in the first 
place, the selfconsistency of the adiabatic approach is guaranteed and leads to localized states in this strong coupling 
limit where the time derivatives of the local deformations can be neglected. As A is varied, these potentials for the 
weak and strong coupling limits are expected to join up smoothly. Depending on the strength of the coupling A, 
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FIG. 1: The variation of V BC (u) (full lines) felt 
by an electron with wave vector k=Q in the pres- 
ence of a static continuous deformation of the in- 
tramolecular distance u for various values of a 
and 7. The dotted line describes the energy of 
the electron with homogeneously deformed oscil- 
lators uncoupled to the electron. The crossing 
points (dots) of these curves indicate a change 
from itinerant to localized behavior. 



this small exercise shows that one can have a situation where the energies of the two configurations are degenerate. 
In this intermediary coupling regime one should expect a coexistence of small localized polarons and weakly bound 
electrons. It has become customary in the theory of small polarons to introduce two dimensionless parameters: 
a = ^e p /hujQ measuring the strength of the interaction and 7 = t/hujo the adiabaticity ratio (a has nothing to do 
with the dimensionless coupling constant ap r used in the theory of large polarons and introduced in section ^J. We 
then can rewrite the potential V sc ({ui}) in a compact form in terms of those parameters and a dimensionless local 
lattice deformation u\ = U{ \J AlLUo/2h. Its variation with a and 7 as a function of u illustrated in fig. ^ 

This picture of the polaron physics, widely used and appreciated until the nineteen sixties (see ref. ? ), although 
intuitively very appealing is qualitatively incorrect - as we shall see in the next section. There we shall show that a 
polaron, although self-trapped, remains delocalized for any dimension. 

A formal and correct treatment of this problem was first initiated by Holstein and collaborators 7 ? and in a very 
elegant and efficient way in terms of a unitary transformation by Firsov and collaborators 7 , given by 

H = e s He- s , S = a^n io .(ai - at) , (30) 

ia 

which transform the electron operators into operators which describe charge carriers rigidly tied to local lattice 
deformations 

= c ia = c ia xr, x ± = e ±^-at) (31) 

and which correspond to localized polaronic self-trapped states 

*i + |0)i = £ e ^ Q2 ^)i- (32) 

n v 

|n)i denotes the n-th excited oscillator states at a molecular sites i and Xj + |0)i signifies an oscillator ground state whose 
equilibrium position is shifted by an amount X/Mujq. The corresponding transformed Hamiltonian is consequently 
given by 

H = J2(D-e p )ctc l<7 -tJ2(4a^X+Xr+H.c.) 

icr i^jcr 

+ ^2(U - 2e p )c+c+c u c iT + tu^o^2(af ai + -) . (33) 

i i 

Such an approach to the polaron problem is particularly useful in the limit of large a where e p > D and for a small 
adiabaticity ratio 7. In the strong coupling limit and anti-adiabaticity (7 < 1), the term X^X^ in the transformed 

Hamiltonian can be averaged over the bare phonon states which, for zero temperature, results in an effective polaron 

2 

hopping integral t* = te~ a . This is justified a posteriori since then t* <C wo, where the local lattice deformations 
can be considered to adapt themselves quasi instantaneously to the slowly in time varying positions of the electrons. 
Such an approximation implies that the number of phonons in the phonon clouds surrounding each charge carrier 
remains largely unchanged during the transfer of a charge carrier from one site to the next, while processes where the 
number of phonons in the clouds change give rise to a polaron damping. Neglecting such damping effects, one obtains 
well defined Bloch states for the polaronic charge carriers (defined by cf a X^) , albeit with a much reduced hopping 
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integral t*, while at the same time the electrons loose practically all of their coherence and thus their quasi-particle 
features. In order to illustrate that, let us consider the Green's function for a single localized polaron, respectively 
for a single localized electron by putting t* = 0. For the polaron retarded Green's function we have 



G p ct (i) = -0(t)(O|X i -[e^ t c+(O)e-^ t c iCT (O)]X+|O) 



which after Fourier transforming becomes 



\ , gi(sp-nftijJo)* 



|(0|X r c+(0)|n)|^ M = -0(t)e^\(0\ct (0)|0)| 



G p et (w) = lim 



1 



5^0 LJ + iS + £r 



(34) 



(35) 



and displays a spectrum which consists of exclusively a coherent contribution. On the contrary, the electron retarded 
Green's function 



G r e f(t) = -e(t)(Q\Xr[e ifit ct(0)e- m c itr (0)}X+\0) 



n=0 



displays a spectrum given by its Fourier transform 



lim ■ 

n=0 



1 



i! lj + id + £ p — rt/iwo 



(36) 



(37) 



The spectral weight of the coherent part is now reduced to exp(— a 2 ), corresponding to the term n — in eq. i|37|l . 
while the major part of the spectrum is made up by the incoherent contributions which track the structure of the 
composite nature of the polaron. 

Generalizing these results by including the itinerancy of the charge carriers leads to very similar results 7 for the 
electron Green's function for a many polaron system, which in this strong coupling anti-adiabatic limit reduces to a 
system of small polarons in band states and where the Green's function for the electrons with wavevectors k is given 
by 



GS t (k,w) = lim 



<5^0 LJ + iS + £ p 



el 



OO 2tj 1 

^ n\ N ^ 

n=l k' 



6->Q 



/(4<) 



1-/(40 



lj + iS - 



+ nfkjQ lj + id — et 



nhuo 



(38) 



and where e£ = e~ a £k- These characteristic spectral properties of small polarons can be tested by photo-emission 
spectroscopy and present crucial tests which permit to distinguishing between different mechanisms leading to heavily 
dressed composite quasi-particles. The Poisson type phonon distribution was particularly well demonstrated early on 
in such experiments on simple molecules such as H2 molecules 7 and carbon rings 7 , on localized polaron states in the 
manganites by neutron spectroscopy 7 and in the cuprate superconductors by infrared absorption measurements 7 . 

The localized, or almost localized, nature of such small polarons is apparent in their electron occupation number 
distribution illustrated in fig. |2 

showing an almost flat distribution covering the whole Brillouin zone, i.e., 

^ = (ckO = (1 ~ e" Q2 K + e" Q V(4) • (39) 

There is a qualitative difference between the weak coupling and the strong coupling limit discussed here. We have 
seen that in the weak coupling limit of the Frohlich Hamiltonian the electrons are accompanied by phonons with which 
they continuously exchange the momentum such as to keep the total momentum of electron plus phonon constant. 
In the strong coupling limit the electrons are surrounded by real lattice deformations which in terms of phonons 
means clouds of phonons and which in the limit of extreme strong coupling, where the electrons can be considered as 
localized on a given site, have a Poisson distribution. The new feature now is that this deformation corresponds to a 
local mass inhomogeneity capable of carrying true momentum rather than the pseudo-momentum as is the case for 
weak coupling. Again, as we shall see below, the momentum between the electron and this deformation is perpetually 
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FIG. 2: Occupation number for a system of 
itinerant small polarons for a>l. The dashed 
and dotted curves indicate the case n = 1 for 
q ~ 1 and a <C 1 respectively (after ref. ' ). 



exchanged as the polaronic charge carrier moves through the lattice, leading to a space-time dependent interaction 
between the charge carrier and the deformation associated to it. This effect will be particularly important in the 
cross-over between the weak coupling adiabatic and the strong coupling anti-adiabatic regime, discussed in the next 
section. 

A key point which occupied the field of polaron theory for several decades was to try to establish if a large 
delocalized polaron changes discontinuously or continuously — although abruptly — into a small self-trapped polaron 
as the electron-lattice coupling is increased beyond a certain critical value. We shall review this issue in the following. 



IV. SELF-TRAPPING 



We shall in this section discuss the meaning of self-trapping. For that purpose let us go back for a moment to the 
semi-classical approach of the polaron problem (discussed in section IIIII above) and the adiabatic lattice potential 
which goes with it. On the basis of such an approach, amounting to treat the polaron problem as a potential problem, 
it is tempting to conclude localization of charge carriers as their coupling to the lattice degrees of freedom increases. 
As we shall see below, the polaron problem can not be treated as a potential problem. Its intrinsic dynamics of the 
local lattice deformations is determinant in correlating the dynamics of the charge and the lattice degrees of freedom, 
thus resulting in itinerant delocalized states of electrons surrounded by clouds of phonons whose density varies and 
increases as we go from the adiabatic weak coupling to the anti-adiabatic strong coupling limit. For that reason the 
lattice degrees of freedom must be treated in a quantized version. 

Let us start with the weak coupling limit for which we can assume a polaron state in the form 



I*)- = 77f E eik ' r, Kk + fM^tm > (40) 



which is the real space version of the state previously discussed, i.e., eq. Q24fl . Diagonalizing the Hamiltonian within 
such a subspace of zero and respectively one phonon present, leads to the following selfconsistent equations for the 
eigenvalues and parameters 



s,-[an, f^ \ (41) 

k' 

1 x - e ik '-( r «- r i) 

fh = ahu -y — . (42) 

k' 

Considering the solution for the ground state (k = 0) of this problem, one notices that its corresponding eigenvalue 
Eq falls below the bottom of the free electron band (Eq < 0). This signals a bound state which is characterized by 
an exponential drop-off of the spatial correlation between the electron and the accompanying phonon, which for a 3D 
system is given by: 

o = avom 1 V2m( _ gn+ ^ o) ^il 
,7lJ 2nh Ir-i-rjl v ; 

(v denoting the unit cell volume) and which indicates that the dynamically (not thermally!) excited phonons accompa- 
nying the electron remain in its immediate vicinity. The exponential drop-off does not depend on any dimensionality 
(in contrast to the semi-classical approach outlined in section IHIjl and gives us a first indication for the intricate 
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correlation between the inherent dynamics of lattice fluctuations and the charge dynamics. Hence, already in this 
weak coupling limit, the polaron problem can not be reduced to an effective (adiabatic) potential problem. 

Let us next turn to the strong coupling limit, where on the basis of such an adiabatic approach one would expect 
self-trapped lowest energy eigenstates of the form 

l*>k = -7=E elk ' ri 4l $ k(K})>, (44) 



where $k({uj}) denotes a function of the ensemble of oscillator coordinates uj situated at sites j in the vicinity of 
site i where the electron is located. Guided by the semi-classical picture of the adiabatic potential, resulting in the 
strong coupling limit, one can try an intuitively appealing Ansatz for the lattice wavefunction in form of a series of 
displaced oscillators around the site where the electron is situated and which, for the ground state, can be assumed 
to be of the form 

|$ k= o({«j})) = e-^J^-^IO) . (45) 

Determining the parameters /y variationally 7 leads to solutions which change discontinuously in the intermediary 
coupling regime. Associated with that is a discontinuous change of the mass of the corresponding quasi-particles 
amounting to a change-over from a practically free band dispersion with a electron mass m c \ to an effective polaron 
mass m p given by m v /m c \ = e a (see fig. 01. 

It has been a matter of dispute for many decades and up to the 1960-ties whether this discontinuity is a real effect or 
is an artifact of the theory. From the experimental side, small polarons where found to be localized and have ever since, 
for the presently available materials, shown a mobility via hopping rather than Bloch like band motion and metallic 
conductivity. The theoretical results obtained for the cross-over regime between large and small polarons on the basis 
of the semi-classical description in terms of an adiabatic lattice potential with two minima of comparable energy (see 
the discussion in section ITTl)l . suggest already that in this regime one should have particularly strong fluctuations of 
the local lattice displacements which could possibly be modeled by some effective lattice wave function in terms of a 
superposition of two sorts of oscillator states: one practically undisplaced and one displaced oscillator state on the site 
housing the electron, as hypothesized early on in a slightly different context 7 . These ideas were subsequently followed 
up by a variety of different approaches such as: exact diagonalization studies 7 on small clusters, upon assuming the 
charge induced deformation to be constraint to the site where the electron sits 7 (when the problem can be solved 
analytically) and Quantum Monte Carlo studies 7 . 

All these approaches converge to a lattice wave function for the ground state which can be approximated by 

I*p({«|})) = ^ aE i^^ _/ii ° j) + jje-^j^-^WQ) . (46) 

This results (after variationally determining the parameters /, g and r/) in a smooth, but nevertheless very abrupt, 
cross-over between the weak and strong coupling features and a change-over from a quasi-free electron band mass to 
a strongly enhanced one (see fig. |3J. The form of eq. (|!l)|> represents for any electron situated on a given site to be 
associated with an oscillator wave function on this site, illustrated in fig. It indicates an essentially undisplaced 
oscillator for weak coupling, a displaced one for strong coupling and a superposition of such two oscillator states for the 
cross-over regime between those two limits. This latter has led to the suggestion of strong retardation effects between 
the dynamics of the electron and that of the local lattice deformation which could possibly result in a dynamically 
disordered system and subsequent localization, a scenario which we shall return to in the lecture "From Cooper-pairs 
to resonating Bipolarons" . 

Let us here consider this point in more detail and investigate this cross-over behavior in terms of the strong 
dynamical correlation which act between the charge and the lattice degrees of freedom. We shall demonstrate that on 
hand of a polaron toy problem, such as a two-site system involving two adjacent diatomic molecules whose individual 
oscillations are uncorrelated with each other and an electron with spin a hopping between those two molecular units. 
The Hamiltonian, eq. I|27|) . for that then reduces to 

H = t(m a + nio) - t{<A<T C 2<r + C ta C l<r) ~ A(niui -I- TI2U2) 

+^-m + ut)+^l + 4)] ■ (47) 



Upon introducing the variables 
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FIG. 3: The effective mass ver- 
sus the electron-lattice coupling u (= 



in our notation) 



after 



ref. ' , for oscillator trial wave functions 
fy 1 and vP such as given by eqs. (1451 . 



FIG. 4: The variation of the oscillator wave function 
on a given site (after ref. ' ) as a function of {u, Wo} 
(= {q 2 (72z)~ 1 ), a 2 } in our notation) for oscillator 
trial wave functions "if 11 such as given by eq. 14611 
and with x — -^(at + en}- 



this Hamiltonian separates into two contributions Hx and Hy depending respectively on lattice coordinates X and 
Y. Y couples exclusively to the total number of charge carriers in the system, ie. (n lcr + n 2a )- It presents an in-phase 
oscillation of the two molecules and has hence no effect on the dynamics of the polaron and will be disregarded for 
the present considerations. The lattice coordinate X, on the contrary, presents an out of phase oscillation of the 
two molecules and couples to the relative charge distribution (ni CT — n 2a ). The Hamiltonian for this subspace which 
describes the intertwined dynamics of the charge and the lattice deformations is 

A A/[ r 

H x = t{n 1<T + n 2a ) - t(c+ a c 2a + c+ c la ) - ^=(n lo - n 2(J )X + y[i 2 + u 2 X 2 \ . (49) 

The term proportional to A in this Hamiltonian clearly indicates how the charge fluctuations between the two molecules 
induces a coupling between the oscillators of those two molecules and which ultimately couples back to the dynamics 
of the charge transfer. This generates, as we shall see, a new dynamical lattice deformation mode which accompanies 
the transfer of charge from one to the other molecule in a slow and smooth continuous fashion. Upon introducing 
the quantification of the out of phase lattice coordinate X = (a + a + )/ y /, 2Mu>o/h in terms of corresponding phonon 
operators a + , a and diagonalizing the above Hamiltonian in a truncated Hilbert space (keeping the number of phonon 
states finite but large) 7 we obtain for the ground state of this system with a single electron with spin a: 

\GS) a = -L [c+M° + {X)) + c + \i>°_(X))] . (50) 

|^(X)) denote oscillator states in real space when the electron is either situated on site 1 or site 
2, similar to the oscillator wave functions illustrated in fig. In particular, in the cross-over 

regime, these oscillator wave-functions reflect the bimodal probability distribution of the equilibrium posi- 
tions of the oscillator where the electron sits. Let us next investigate the evolution in time of the 
charge transfer together with that of the inter-molecular deformation, given by the correlation functions: 

X nn (r) = 9(r)([n 1<T (r) - n 2r7 (r)][n lr7 (0) - n 2CT (0)]) (51) 



X xx =6(t)(X(t)X(0)). (52) 

The results of this are reproduced in fig. [S] for a fixed value of a (which in the present notation corresponds to 
a = y2 • 1.2 = 1.70) and where we cover the cross-over from the strong coupling anti-adiabatic limit (7 = 0.1) to the 
strong coupling adiabatic limit (7 = 2.0) upon increasing the adiabaticity ratio. 
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FIG. 5: The evolution in time (in units of uIq 1 ) 
of the charge and intra-molecular distance trans- 
fer between two adjacent molecular units for a 
two-site one-electron system (after ref. ' ) in the 
anti-adiabatic (7 = 0.1), the adiabatic (7 = 2.0) 
and in the cross-over regime (7 = 1.1) and for 
(in the present notation) a coupling strength 
a = 1.2V2 ~ 1.7. 



The anti-adiabatic limit is characterized by a smooth and slowly in time varying transfer of charge from one 
molecule to the next. This charge transfer then slaves the inter-molecular deformation X by subjecting it to a slowly 
in time varying driving force — -^[ni^r) — n2 ff (r)]X(r). This in turn leads to a slowly in time varying sinusoidal 
intermolecular deformation, onto which are superposed the local intrinsic oscillations of the individual molecules with 
a frequency of the order of u)q. This is the exact opposite of what happens in the extreme other limit, when we 
approach the adiabatic situation. 

In this adiabatic regime, considering here an adiabaticity ratio 7 = 2.0, the slowly in time varying quantity is 
now the inter- molecular deformation (showing almost no fluctuations of the individual molecules with frequency wo), 
which slaves the intermolecular charge transfer characterized by [ni CT (r) — rc.2o-(r)] • Averaged over several periods of 
oscillations of frequency ~ 1/t, this charge transfer adiabatically follows again the slowly in time varying potential 
-^[n la (r)-n 2a {T)]X{T). 

These results show that in these two extreme limits of adiabaticity and anti-adiabaticity the slowly in time varying 
component of the charge transfer processes and of the inter-molecular deformation dynamics are completely locked 
together. This is no longer the case in the cross-over regime between these two extreme limits, here shown in fig. 
for the adiabaticity ratio 1.1. One can clearly distinguish phase-slips occurring between the correlated motion of the 
charge and the deformation transfer, which, in a large system where many body effects become important, could 
possibly lead to a dynamically induced localization of the charge carriers. 



V. SOME EARLY EXPERIMENTAL INSIGHTS 



The local character of small polaronic charge carriers requires specific experimental probes able to track the polaron 
induced lattice as well as charge excitations on a short spatial (~ 5 A) as well as time (10~ 13 — 10~ 15 sec) scale. 
Generally such probes are quite adequate for insulating polaronic systems with either charge ordered bipolarons 
(Ti40y), bipolaronic Mott type insulators (WO3-U the manganites and nickelates) and low density polaron systems 
arising from photo-induced or chemical doping and situated on the verge of a quantum phase insulator-superconductor 
transition (high T c cuprates). Concerning the cuprates, manganites and nickelates we refer the reader to the lectures 
by T. Egami and N. L. Saini in this volume, dealing with pulsed inelastic neutron scattering techniques and EXAFS 
as well as XANES. Here we shall restrict ourselves to the discussion of more classical probes such as optical absorption 
which can select specific local polaron-sensitive lattice modes. 

As we have seen in the discussion presented above, one of the major issue here is to track the disintegration and 
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FIG. 6: Time evolution of the number of W 
sites (at different concentrations created with 
varying pump laser intensities breaking the 
bipolaronic W 5+ -W 5+ units) after after shut- 
ting off the pump-laser (after ref. ' ). The full 
lines represent the theoretical bimolecular re- 
laxation behavior. The inset is a schematic 
picture for the bimolecular relaxation process. 
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FIG. 7: The photoinduced infrared absorp- 
tion coming from the laser-stimulated, re- 
spectively laser-suppressed vibrational lattice 
modes (low frequency part) and the photo- 
modulation response for electrical ingap exci- 
tations (high frequency part) in YBa2Cu30e.2 
(after ref.' ). 



reconstruction of a polaron during its transfer from one site to the next. This process implies a gradual stripping 
off of the electron's phonon cloud upon leaving a given site and a subsequent rebuilding of this phonon cloud on the 
new site where the electron eventually ends up. Another important issue in the polaron problem is connected to 
the physics of Many-Polaron system, caused by the polaron induced local attraction which can bind two polarons 
on a given effective site into a bipolaron. We shall now discuss specific experiments which can illustrate these two 
characteristic polaron features. 

Shortly after Anderson's suggestion 7 of localized bipolarons in amorphous chalcogenide glasses, exhibiting a natural 
diamagnctism resulting from covalent bonding in locally deformed lattice structures, a variety of systems were found 
where bipolarons existed in more or less dense situations. Examples for dense bipolaronic systems, exhibiting spatially 
symmetry broken states related to bipolaron ordering are: Ti^^V^Oy 7 and Naa;V205 ? . An example for a dilute 
bipolaron systems is W03_ x ? . In those systems the bipolarons form on adjacent cations, such as Ti 3+ -Ti 3+ , V 4+ -V 4+ 
and W 5+ -W 5+ bonds in strongly deformed octahedral ligand environments. Those bipolaronic units are imbeded in a 
corresponding background of Ti 4+ -Ti 4+ , V 5+ -V 5+ and W 6+ -W 6+ molecular units, which together with the bipolaronic 
entities, constitute those crystalline materials. 



A. Bipolaron dissociation and recombination in W03- T 

Bipolarons can be dissociated with light which, in the case of WOs-^, leads to the creation of isolated W 5+ cation 
sites in vibrationally excited states of the ligand environments and whose concentration can be tracked by electron 
spin resonance (ESR) signals coming from those W 5+ sites. Optical absorption measurements 7 in non-illuminated 
samples show a peak in the spectrum centered around 0.7 eV coming from pre-existing isolated W 5+ cations sites. 
When the crystal is illuminated with a broad spectral band centered around a suitable frequency (1.1 eV in this 
case) the optical absorption coming from the W 5+ cations sites increases while, concomitantly, an intrinsic shoulder 
of the absorption band at around 1 eV — attributed to the absorption coming from the intrinsic bipolaronic W 5+ -W 5+ 
bonds — decreases correspondingly. After the illumination is shut off, the single-polaron vibrationally exited W 5+ 
units relax to bipolaronic ones, as can be tracked by ESR measurements as a function of time (see fig-EJ. 

The time evolution of this relaxation follows a bimolecular recombination process controlled by a rate equation 
given by 

= -B(t)n 2 w , (53) 

where nw denotes the density of the W 5+ ions and B(t) is a function describing diffusion limited reactions. This 
relaxation behavior is clearly distinct from that of mono-molecular process which describe relaxation processes from 
vibrationally excited single polaron units to their ground state. Such time resolved studies permit to investigate the 
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rebuilding of the phonon cloud for a bipolaron after its separation into two separate polarons, much as would be 
expected for processes of bipolarons hopping between neighboring sites. 

B. Photo-induced polarionic charge carriers in high T c cuprates 

A somewhat related in spirit technique was used in the study of polaronic features in the high T c cuprates 7 . 
It was based on examining the variation in the optical absorption spectrum upon doping those materials in the 
semiconducting phase with electron-hole pairs using so-called photo-induced doping. The idea was to use laser 
illumination with photon energies bigger than the semiconductor gap in the pump laser beam. In the specific case of 
the high T c cuprates this meant a photon energy of typically 2 to 2.5 eV and which created an estimated density of 
charge carriers of about 10 19 per cm 3 . The injected photoinduced charge carriers modify locally the lattice symmetry 
and, by doing so, lead to the activation of corresponding local phonon modes (the 434 cm' 1 out of plane mode and 
the 500 cm -1 axial [CU-O4] stretching mode). These modes can be made evident as the steady state response of the 
system to the photo illumination in the absorption spectrum tested by a probe infrared beam (see fig.0. In such an 
experimental set up, the 590 cm -1 mode with the tetragonal symmetry (being associated with the undistorted lattice 
which characterize the non-illuminated materials) loose in intensity (bleaching effect) in the absorption spectrum 
after the samples have been illuminated. Examining the response of the system at higher frequencies, covering the 
energy range of the semiconducting gap and higher, by photo-modulation techniques (which test the photo-induced 
charge carriers by their absorption of light in a narrow frequency window) visualizes the excitations with a life time 
which is inversely proportional to this frequency. The corresponding absorption spectra show an activation of charge 
excitations over a broad background, covering the energy region of the semiconducting gap, and a bleaching for 
energies above that frequency (see fig. EJ. These results imply a shift of spectral weight upon illumination from the 
states above the gap in the non-illuminated samples into in-gap states. The intensity of the absorption of the activated 
phonon modes as well as of the activated electronic in-gap excitations turns out to scale like the square root of the 
intensity of the laser pulse illuminating the sample. This suggests that (in analogy with studies on WO^^ x discussed 
above) those excitations relax via bimolecular recombination processes. Since the studies on photo-induced charge 
carriers have features similar to those in low doped systems, obtained by chemical substitution, it has been tempting 
to conclude that the charge carriers in the low doped superconducting samples have resonant bipolaronic features of 
a finite lifetime rather than corresponding to well defined stable bipolaronic entities. Such entities can nevertheless 
condense into a superconducting state and are controlled by phase rather than amplitude fluctuations as discussed in 
my lecture: "From Cooper pairs to resonating bipolarons" in this volume. 

C. Bipolaronic charge ordering in Ti4C>7 

A particularly interesting and physically very rich example of a dense bipolaronic system is Ti407. It exhibits a low 
temperature phase (T < 140 K) where bipolarons are in a symmetry broken ordered state, consisting of diamagnetic 
T 3+ — T 3+ bipolaronic diatomic pairs with sensibly reduced intra-molecular distances as compared to T 4+ — T 4+ 
pairs with which they alternate in quasi 2D slab like structures. Upon increasing the temperature, there is a small 
interval, [140 K < T < 150 K], were these bipolaronic electron pairs are dynamically disordered and eventually break 
up into individual electrons and a metallic phase for T > 150 K. The experimental measurements of this material 
involved x-ray diffraction, resistivity, susceptibility and specific heat measurements as well as electron paramagnetic 
resonance (EPR) studies, summarized in ref. ? . The two low temperature phases in the regimes T < 140 K and 140 
K < T < 150 K are semiconducting with a conductivity characterized by comparable activation energies of the order 
of 0.16 eV. The phase transition between those two semiconducting phases can be attributed to an order-disorder 
phase transition where the bipolarons are essentially the same as in the low temperature ordered phase, as evident 
from the absence of any change in the magnetic susceptibility as well as intra-molecular distance of the bipolarons 
when going through this phase transition. The transition to the high temperature phase at T = 150 K is characterized 
by a breaking up of the bipolarons into itinerant electrons leading to a metal with an enhanced Pauli susceptibility 
and a substantial decrease in the unit-cell volume (see fig.|SJ). Both transitions are first order and the entropy change 
in the high temperature transition is due, in roughly equal amounts, to electronic and lattice contributions. The 
most important feature of this system is however that the disordered bipolaronic phase is a dynamical rather than 
static disorder as evidenced from EPR experiments which show a vanishing of the EPR line upon entering this phase 
from the low temperature ordered phase. It was on the basis of these experimental findings that the possibility of a 
condensation of bipolarons was initially proposed 7 . This broke with a traditionally severely guarded doctrine and 
stimulated to look for superconducting materials which (i) were oxides, (ii) have reduced dimensionality, (iii) are close 
to insulating parent compounds and (iv) are generally poor rather than good metals in the normal phase. 
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FIG. 8: The electrical conductivity for Ti 4 07 
(after ref. ' ) showing with increasing T a tran- 
sition from a low temperature activated hopping 
regime to a different activated hopping regime at 
T ~ 140 K followed by a second transition from 
that to a metallic phase at T ~ 150 K. The mag- 
netic susceptibility, shown in the inset, does not 
change at the low temperature transition 



lOoo/nK" 1 ] 



On the basis of such a scenario and the Hostein model one could expect for such systems a superfluid state of 
bipolarons, albeit with a very small critical temperature, since being inversely proportional to the bipolaron mass 
which is typically several orders of magnitude bigger than the electron band mass. The real difficulty to observe this 
type of superfluidity in crystalline materials might however be related to the fact that the standard polaron models, 
generally based on harmonic lattice potentials, totally neglect any relaxation processes. High T c cuprates are clearly 
not candidates for this extreme case of Bipolaronic Superconductivity but are likely to contain localized bipolarons 
as resonant states inside the Fermi sea of itinerant electrons which could result in a superconducting state controlled 
by phase rather than amplitude correlations. 



VI. SUMMARY 



In this introductory lecture I briefly reviewed various kinds of electron-lattice couplings, characterizing different 
classes of materials and which give rise to two distinct categories of polarons: large Frohlich and small Holstein 
polarons. The importance of treating the phonons as quantum rather than classical variables became evident in 
connection with the question of itinerancy of the polaronic charge carriers. The dynamics of the polaron motion, 
exemplified in real time, shows a highly non-linear physics involving the dynamics of coupled charge and the lattice 
fluctuations which mutually drive each other. The present discussion was restricted to the single polaron respectively 
bipolaron problem and to limiting cases (such as strong coupling anti-adiabatic limit) where the Many Polaron problem 
can be decomposed into a band of single polaron states with different wave vectors. The fundamental questions of 
the polaron problem which pose themselves today, are evidently beyond the topics touched upon in this preliminary 
discussion. These are questions which concern the cross-over between the adiabatic and anti-adibatic regimes, the 
polaron induced residual interactions in a Many Polaron system and its dependence on the range of electron-lattice 
coupling as well as on the density of charge carriers, which can result in possible transitions between insulating and 
metallic behavior of polarons. At this stage, in order to tackle this kind of problems we have to resort to highly 
sophisticates numerical techniques, which will be presented in this school. The hope is that eventually this will give 
us some insight into this complex Many Body problem so that sooner or later we can formulate this polaron physics 
in a way where analytical approaches can capture its main qualitative features. 
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